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Allowed regions of core fraction versus partial coherence are obtained from analysis of 
both two- and three-pion Bose-Einstein correlations measured at PHENIX in ^/snn = 200 
GeV Au+Au collisions. Two-pion Bose-Einstein correlation functions for different average 
transverse momenta are used to analyze chiral symmetry restoration. 
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1. Introduction 

q . Correlation functions are important to see the collective properties of particles and 

space-time structure of the emitting source, eg. the observed size of a system can be 
measured by two particle Bose-Einstein correlations [ . The m t dependent strength 



of two-pion correlations provides information about Ua(1) symmetry restoration in the 
source [0EHUE1- From a consistent analysis of two- and three-particle correlations one 
can establish an experimental measure of thermalization and coherence in the source [E]. 

2. Data Analysis 



We studied 70M v /5^ = 200 GeV RHIC Au+Au collisions recorded at PHENIX during 
Run 4. Charged tracks were detected by the Drift Chamber and Pad Chambers. The 
Time Of Flight detector and the Electromagnetic Calorimeters were used to identify 7r + 
particles, and we had 200M 7r + 's, 900M tt + pairs and more than 4G 7r + triplets in this 
analysis. We measured from this data sample the two-pion correlation function C 2 as a 
function of qw = \J {k\ — k 2 ) 2 and the three-pion correlation function C 3 as a function of 
q inv of the three pairs in the triplet, qi 2 , q 2 3, <?3i- 

Our method of Coulomb- correct ion is based on the solution of the two-body Coulomb- 
problem. A symmetrized asymptotically correct three-body Coulomb wave-function is 
built up from the two-body solution. The source function is then integrated with the 
absolute square of this three-body Coulomb-wave function as the density function, and 
the same is done for a symmetrized three-body plain-wave function. The ratio of the two 
gives the Coulomb- correct ion factor [Cj. For the final Coulomb-correction we utilized a 
core/halo type of picture of the source and iterated A and R in a self-consistent manner 
until convergence. 

The fits were done using three different shapes, Gauss, Levy and Edgeworth, described 
in ref. [|H]. We fitted only those points with q inv > 0.02 GeV/c due to a non-BEC and 
non-Coulomb structure seen at low q inv (see the gray lines on fig. |TJ). 



*For the full list of PHENIX authors and acknowledgments, see appendix 'Collaborations' of this volume. 
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Figure 1. Raw- and Coulomb-corrected correlation functions for 0.2 < p t < 2.0 GeV/c 
and — 92% centrality with Levy fits. Points with q inv < 0.02 GeV/c were excluded from 
the fits. 



3. Results 

First, we fitted two- and three-particle correlation functions for 0.2 < p t < 2.0 GeV/c 
and 0—92% centrality. Here only the Levy fits had acceptable confidence levels. These and 
the Coulomb-corrected correlation functions are shown on fig. ^ The three-dimensional 
C 3 was projected onto the q 12 = q 23 = Q31 line. 

Let us consider the ratios f c = N c (p)/Ni(p) (fractional core) and p c = N^{p)/N c {p) 
(partial coherence), and take the equations (see eqs. 57 and 58 in ref. [|S]): 

C 2 ( Pl ^ P2 ) = l + f c [l-p 2 c } (1) 
C 3 ( Pl ~p 2 ~p 3 ) = i + 3/ c 2 [l-^] + 2 / c 3 [(l-p c ) 3 + 3p c (l-p c ) 2 ] (2) 

From this we calculated the allowed f c vs p c region, using the 3a contours from eqs. 
121 see fig. 121 The size of the region obtained from C 3 is much larger due to the smaller 
number of triplets compared to pairs. 

Compared to former NA44 analysis (see refs. [El El El), the resulting allowed region is 
at a little bit higher f c values, because our f c was obtained from Levy fits and not from 
Gauss, as confidence level of latter was unacceptable. 

In the second part of the analysis, we have analyzed two-pion correlation functions in 
10 pt bins from 0.2 GeV/c to 0.5 GeV/c with — 92% centrality, using the invariant 
momentum difference variable and the following shapes (in he = 1 units): 

Gauss: C 2 (q) = 1 + Aexp(- \Rq\ 2 ) (3) 
Levy: C 2 (q) = 1 + Aexp(- \Rq\ a ) (4) 

Edgeworth: C 2 {q) = 1 + Aexp(- |i2g| 2 )(l + ^H 3 (y/2 \Rq\), (5) 

where H 3 is the third Hermite polynomial, and k 3 it's coefficient (the Edgeworth form is a 
model-independent expansion, of which we keep only the first term, to measure deviation 
from the Gauss form). See more details in ref. [Illj. 

The parameters as a function of average transverse mass (A(m t ), R(m t ), a(m t ) and 
K 3( m t)) for Gauss, Levy and Edgeworth fits are shown on fig. 01 The a parameter of the 
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Figure 2. Allowed region of partial coherence versus fractional core extracted from Levy 
fits. Contours obtained from C3 are much larger due to smaller statistics of 7r + triplets 
compared to n + pairs. 

Levy distribution controls the long-range source in a similar way as the A parameter, so 
fitting both makes the fit underconstrained. Hence we did fits also with fixed a, and for 
comparison, also with fixed k 3 . These latter fits show the same \{m t ) behavior as in the 
Gaussian case. 

The m t dependence of A can be used to extract information on the mass-reduction of 
the rf meson, a signal of the Ua(1) symmetry restoration in hot and dense matter (see 
ref. [ 2J). A comparison of the measurements of fig. H3 with model calculations of ref. 
using Fritiof results for the composition of the long-lived resonances and a variation of 
the rf mass is presented in fig. HJ 

4. Conclusion 

From simultaneous analysis of two- and three-pion correlations we obtained an allowed 
region of partial coherence versus core fraction. Nonzero partial coherence is allowed when 
having larger core size, at total incoherence (p c = 0) the core fraction is slightly higher 
than previous analysis using Gaussian fits (see refs. [0EIIIDIII21)- This is due to the fact 
that we used Levy fits, the only investigated one that gave acceptable confidence levels. 

We also analyzed the m t dependence of the A parameter obtained from various fits. 
Gauss fit results agree with former PHENIX measurements (see ref. Regarding 
Ua(1) symmetry restoration, we conclude that at present, results are critically depen- 
dent on our understanding of statistical and systematic errors, and additional analysis is 
required to make a definitive statement. 
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Figure 3. A(mt), R(mt), ct(mt) and K^{mt) from Gauss, Levy and Edgeworth fits 
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Figure 4. Measured \{m t ) compared to calculations using the model of ref. [Ej with various 
rf mass values, and the Fritiof results for the composition of the long-lived resonances. 
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